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Abstract 

The cosmological bound of the stochastic background of gravitational 
waves is analyzed with the auxilium of the WMAP data, differently from 
lots of works in literature, where the old COBE data were used. From 
our analysis, it will result that the WMAP bounds on the energy spec- 
trum and on the characteristic amplitude of the stochastic background of 
gravitational waves are greater than the COBE ones, but they are also far 
below frequencies of the earth-based antennas band. At the end of this 
letter a lower bound for the integration time of a potential detection with 
advanced LIGO is released and compared with the previous one arising 
from the old COBE data. Even if the new lower bound is minor than the 
previous one, it results very long, thus for a possible detection we hope 
in the LISA interferometer and in a further growth in the sensitivity of 
advanced projects. 

The design and construction of a number of sensitive detectors for gravita- 
tional waves (GWs) is underway today. There are some interferometers like 
the VIRGO detector, being built in Cascina near Pisa, Italy, by a joint Italian- 
French collaboration, the GEO 600 detector being built in Hanover, Germany 
by a joint Anglo-Germany collaboration, the two LIGO detectors being built 
in the United States (one in Hanford, Washington and the other in Livingston, 
Louisiana) by a joint Caltech-Mit collaboration, and the TAMA 300 detector, 
being built near Tokyo, Japan. There are many bar detectors currently in op- 
eration too, and several interferometers and bars are in a phase of planning and 
proposal stages (for the current status of GWs experiments see [UH]). 
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The results of these detectors will have a fundamental impact on astrophysics 
and gravitation physics. There will be several experimental data to be analyzed, 
and theorists will be forced to interact with lots of experiments and data analysts 
to extract the physics from the data stream. 

In this letter the stochastic background of GWs (SBGWs) 02IH02, which is, 
in principle, a possible target of these experiments, is analyzed with the auxilium 
of the WMAP data [6l [7] . We emphasize that, in general, in previous works in 
literature about the SBGWs, old COBE data were used (see 02 H 02 03 M for 
example) . 

From our analysis, it will result that the WMAP bounds on the energy 
spectrum and on the characteristic amplitude of the SBGWs are greater than 
the COBE ones, but they are also far below frequencies of the earth-based 
antennas band. At the end of this letter a lower bound for the integration time 
of a potential detection with advanced LIGO is released and compared with the 
previous one arising from the old COBE data. Even if the new lower bound is 
minor than the previous one, it results very long, thus for a possible detection 
we hope in the LISA interferometer and in a further growth in the sensitivity 
of advanced projects. 

The strongest constraint on the spectrum of the relic SBGWs in the fre- 
quency range of ground based antennas like bars and interferometers, which is 
the range 10Hz < f < lO i Hz, comes from the high isotropy observed in the 
Cosmic Background Radiation (CBR). 

The fluctuation AT of the temperature of CBR from its mean value To = 
2.728 K varies from point to point in the sky [6j [T] , and, since the sky can be 
considered the surface of a sphere, the fitting of AT is performed in terms of a 
Laplace series of spherical harmonics 

AT °° 1 

— (0) =EE ai m Y lm (Cl), (1) 

1=1 m=-l 

and the fluctations are assumed statistically independent (< a; m >= 0, 
< ai m a^, m , >= CiSu'Smm'). In eq. ([1]) ft denotes a point on the 2-sphere while 
the ai m are the multipole moments of the CBR. For details about the defini- 
tion of statistically independent fluctations in the context of the temperature 
fluctations of CBR see [10] . 

The WMAP data [6l [7] permit a more precise determination of the rms 
quadrupole moment of the fluctations than the COBE data 

Q rms ^T(J2 l - S f^) i = 8 ± ^K, (2) 

m=-2 

while in the COBE data we had [3J [T3] 

Q rms = U.Stl^K- (3) 

A connection between the fluctuation of the temperature of the CBR and the 
SBGWs derives from the Sachs- Wolfe effect [HOI]. Sachs and Wolfe showed that 
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variations in the density of cosmological fluid and GWs perturbations result 
in the fluctuation of the temperature of the CBR, even if the surface of last 
scattering had a perfectly uniform temperature [H] • In particular the fluctuation 
of the temperature (at the lowest order) in a particular point of the space is 



The integral in eq. (|4|) is taken over a path of null geodesic which leaves the 
current spacetime point heading off in the direction defined by the unit vector 
O and going back to the surface of last scattering of the CBR. 

Here A is a particular choice of the affine parameter along the null geodesic. 
By using conformal coordinates, we have for the metric perturbation 



and r in cq. ((31) is a radial spatial coordinate which goes outwards from the 
current spacetime point. The effect of a long wavelenght GW is to shift the 
temperature distribution of CBR from perfect isotropy. Because the fluctations 
are very small (< AT/Tq >< 5 * 10 -5 0(7]), the perturbations caused by the 
relic SBGWs cannot be too large. 

The WMAP results give rather tigh constraints on the spectrum of the 
SBGWs at very long wavclenghts. In [3] we find a constraint on Q gw (f) which 
derives from the COBE observational limits, given by 



Now the same constraint will be obtained from the WMAP data. Because 
of its specific polarization properties, relic SBGWs should generate particular 
polarization pattern of CBR anisotropies, but the detection of CBR polariza- 
tions is not fulfiled today [T2]. Thus an indirect method will be used. We 
know that relic GWs have very long wavelenghts of Hubble radius size, so the 
CBR power spectrum from them should manifest fast decrease at smaller scales 
(hight multipole moments) . But we also know that scalar modes produce a rich 
CBR power spectrum at large multipole moments (series of acoustic peaks, ref. 
[HI EJ) ■ Then the properties of tensor modes of cosmological perturbations of 
spacetime metric can be extract from observational data using angular CBR 
power spectrum combined with large scale structure of the Universe. One can 
see (fig. 1 ) that in the range 2 < / < 30 (the same used in [3J, but with the 
old COBE data [13]) scalar and tensor contributions are comparable. From 
[51 [7J the WMAP data give for the tensor /scala ratio r the constraint r < 0.9. 
(r < 0.5 in the COBE data, ref. [13]); Novosyadly and Apunevych obtained 
Qscaiar(Ho) < 2.7 * 10~ 9 [12] . Thus, if one remembers that, at order of Hubble 
radius, the tensorial spectral index is —4 < n t < —2 , it results 





5 g a b = R 2 (r])h ab , 



(5) 



fWi>?00 < 7 * 10- n (-^) 2 for H < f < 30H . 



(0 



n g w(f)h 2 100 < 1.6 * HT 9 (-^) 2 for Hq < / < 30ff , 



(7) 
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Spherical Harmonics I 



Figure 1: The tensor to scalar ratio 



which is greater than the COBE data result of eq. ((6]). 

We emphasize that the limit of eq. ([7]) is not a limit for any GWs, but only 
for relic ones of cosmological origin, which were present at the era of the CBR 
decoupling. Also, the same limit only applies over very long wavelenghts (i.e. 
very low frequencies) and it is far below frequencies of the Virgo - LIGO band. 

The primordial production of the relic SBGWs has been analyzed in [3J IH 
[51 U3J US] , where it has been shown that in the range 10~ 15 Hz < f < 10 10 Hz 
the spectrum is flat and proportional to the ratio 

Pd> ~ 10 -12 (8) 



PPlanck 

WMAP observations put strongly severe restrictions on the spectrum, as 
we discussed above. In fig. 2 the spectrum £l gw is mapped: the amplitude 
(determined by the ratio ^ ds - ) has been chosen to be as large as possible, 
consistent with the WMAP constraint ([7]). Nevertheless, because the spectrum 
falls off oc f~ 2 at low frequencies [3JIH[5], this means that today, at Virgo and 
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LISA frequencies, indicated in fig. 2, 



fW/)^ioo O * 10- 13 , (9) 
while using the COBE data it was 

*W/)kioo < 8 * 10" 14 (refs. [31II3]). 

It is intersting to calculate the correspondent strain at rj 100Hz, where 
interferometers like Virgo and LIGO have a maximum in sensitivity. The well 
known equation for the characteristic amplitude [3] [5] can be used: 

h c (f) ~ 1.26 * 10- 18 (^)y^V(7), (10) 

obtcning 

/i c (100Hz) < 1.7* 10" 26 . (11) 

Then, because for ground-based interferometers a sensitivity of the order of 
10~ 22 is expected at rj IQQHz, four order of magnitude have to be gained in 
the signal to noise ratio [TH1 H7]- Let us analyze smaller frequencies too. The 
sensitivity of the Virgo interferometer is of the order of 10~ 21 at « 10Hz [16] 
and in that case it is 



h c (100Hz) < 1.7* 10~ 25 . (12) 

For a better understanding of the difficulties on the detection of the SBGWs 
a lower bound for the integration time of a potential detection with advanced 
LIGO is released. For a cross-correlation between two interferometers the signal 
to noise ratio (SNR) increases as [3] [4] [5j [18] 

™^lif*Sr (I3) 

where Pi{\f\) is the one-sided power spectral density of the i detector [18] 
and j(f) the well known overlap-reduction function [TU [19]. Assuming two 
coincident coaligned detectors (7(/) = 1) with a noise of order 10 _48 /iJz (i.e. 
a typical value for the advanced LIGO sensitivity [20]) one gets (SNR) ~ 1 for 
~ 3* 10 5 years using our result VL gw (f)h\ w ~ 9* 10~ 13 while it is (SNR) ~ 1 for 
~ 3 * I0 7 years using previous COBE result flg W (f)hf Qa ^8* 10~ 14 . Since the 
overlap reduction function degrades the SNR, these results can be considered 
a solid upper limit for the advanced LIGO configuration for the two different 
values of the spectrum. 

The sensitivity of the LISA interferometer will be of the order of 10 -22 at 
10~ 3 w Hz [21] and in that case it is 
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h c (100Hz) < 1.7* 10 -21 . 



(14) 



Then a stochastic background of relic gravitational waves could be in prin- 
ciple detected by the LISA interferometer. We also hope in a further growth in 
the sensitivity of advanced projects. 

We emphasize that the assumption that all the tensorial perturbation in the 
Universe are due to a SBGWs is quit strong, but our results (|9]), ([TTjl . (fT2|) and 
(fl"4"|) can be considered like upper bounds. 

Rcasuming in this letter the SBGWs has been analyzed with the auxilium of 
the WMAP data, while previous works in literature, used the old COBE data, 
seeing that the predicted signal for these relic GWs is very weak. From our 
analysis it resulted that the WMAP bound on the energy spectrum and on the 
characteristic amplitude of the SBGWs are greater than the COBE ones, but 
they are also far below frequencies of the earth-based antennas band. In fact the 
integration time of a potential detection with advanced interferometers is very 
long, thus, for a possible detection we have to hope in a further growth in the 
sensitivity of advanced ground based projects and in the LISA interferometer. 
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Figure 2: The spectrum of relic gravitons in inflationary models is flat over a 
wide range of frequencies. The horizontal axis is log 10 of frequency, in Hz. The 
vertical axis is log 10 VL gr . The inflationary spectrum rises quickly at low frequen- 
cies (wave which rentered in the Hubble sphere after the Universe became mat- 
ter dominated) and falls off above the (appropriately redshifted) frequency scale 
fmax associated with the fastest characteristic time of the phase transition at 
the end of inflation. The amplitude of the flat region depends only on the energy 
density during the inflationary stage; we have chosen the largest amplitude con- 
sistent with the WMAP constrain discussed earlier: il gr (f)hi 00 < 1.6* 10 -9 at 
10~ 18 iiz. This means that at Virgo and Lisa frequencies, ^l gr {f)h\ 00 < 9*10~ 13 
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